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Abstract. We introduce a new class of stochastic processes which are 
defined on a random set in W^. These processes can be seen as a link be- 
tween random fields and marked point processes. Unlike for random fields, 
the mark covariance function need in general not be positive definite. This 
implies that in many situations the use of simple geostatistical methods ap- 
pears to be questionable. Surprisingly, for a special class of processes based 
on Gaussian random fields, we do have positive definiteness for the corre- 
sponding mark covariance function and mark correlation function. 

Classification. Primary: 60G60, 60G55; secondary: 60G15, 60D05 

Keywords, random field, random set, marked point process, mark cor- 
relation function, mark covariance function 

1 Introduction 

Quantities measured in are mostly modelled as so-called regionalized 
variables, i. e., one usually assumes that these quantities can, in principle, be 
measured everywhere in and that the choice of sampling points does not 
depend on the values of these quantities. Based on this assumption, several 
geostatistical methods like variogram analysis or kriging can be applied, see 
[6(1 . However, there are two types of situations where this assumption does 
not hold and hence, uncritical use of geostatistical methods might cause 
incorrect or meaningless results. 
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The first type of problems is caused by the investigators themselves by 
some kind of preferential sampling For instance, this happens when data 
are sampled only at places where high values of the variable of interest are 
expected. The second type of problems is intrinsic to the investigated object 
itself. An obvious situation is the investigation of individuals, e. g. trees in a 
forest, where interactions among individuals are present. In this particular 
situation the theory of marked point processes provides a formal framework 
for data analysis 13, 1^ . 



Here, we like to draw the attention to some further, deceptive situations 



where implicit conditioning has been mostly ignored in literature |l3l.ll6l.l29l|. 
For instance, the investigation of pesticides in soil is restricted to cropland 
and the height of forest litter is restricted to silvicultural areas. In both cases, 
a preselection cannot be excluded since environmental conditions directly 
influence the kind of land use. A further, simple example has motivated this 
work and appears when the altitude is predicted by geostatistical methods 
based on measurements that are taken above sealevel only. 

Such kind of conditioning might be considered as minor, but can cause 
major effects, nonetheless. We advise caution because of the following facts: 

1. Any characteristic, such as the covariance function or the variogram, 
has to be understood as a conditional quantity given measurements can be 
taken at certain locations. 

2. In general, neither the covariance function is positive definite nor the 
variogram is conditionally negative definite. 

Since Gaussian random fields are rather popular, a bigger part of this 
paper deals with the following model: the sealevel is at and the altitude is 
given by some (smooth) stationary Gaussian random field Z with mean —t 
and variance 1. Then we face the following oddities when inference is based 
on measurements above sealevel only: 

1. The theoretical variogram is not conditionally negative definite, in 
general. 

2. A naive definition of the covariance function C{x,y) by 

C(x, y) = E[Z{x)Z{y) \ Z{x) > 0, Z{y) > 0] - rrl^ 

leads in general to a function which is not positive definite, for any m G M. 

3. A more suitable definition of the covariance function for the altitude 
above sealevel as the conditional covariance given that Z{x) > and Z{y) > 
leads to a function which is never differentiable. 

4. If t = 0, the conditional covariance function is positive definite. 
Though, no random field exists that is independent of the sampling loca- 
tions and that can model the altiditude above sealevel. 
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Before discussing the above set up in detail, we will introduce a theo- 
retical framework so that both a meaningful definition of second-order char- 
acteristics is possible and usual random fields as well as marked point pro- 
cesses are included as particular cases. For this reason we extend the no- 
tion of a random upper semi-continuous (u.s.c.) function (taking values in 
M = [— cxD,oo]) on R'^ such that the domain is a random subset of M"^. To 
this end we make use of Matheron's fv^ idea and consider the hypograph 

Af = {{x, t) £ X xR: t< /(x)}, X OR'^ 

of a function f : X ^ R. In fact, ^/ is closed if and only if / is u. s. c. on 
closed X, and the mapping f Af is a bijection. 

The paper is organized as follows. In Section [2] we formally introduce 
the notion of a random marked closed set and discuss some examples. In 
Section [3] we generalise the definition of several characteristics for random 
fields to random marked sets. We show that, in general, they do not share 
the same definiteness properties as their random field analogues. In Section 
m we study Gaussian random fields Z{x) given that Z{x) exceeds a certain 
threshold i € M. In Section [5] some results on the differentiablity of the mark 
covariance function of random marked sets are given. In Section[6l we collect 
the proofs of the statements of the preceding sections. 

2 Random marked closed sets 

Denote by M = M U {— oo, -|-oo} the extended real line. Let 

= {{XJ) : X QR'^ is closed, / : X ^ 1 is u. s. c.}. 

^usc is isomorphic to the system Ud of all closed sets ^ C M'* x M which 
satisfy 

Vx G M'^ Vt G M : (x, t) G yl ^ {x} x [-oo, t] C A (1) 
by the bijection 

T : ^usc Ud 
{X, f) ^ {(x, t)eXxR:t< fix)}, (X, /) G ^usc- 

The subsequent proposition follows immediately from the fact that the space 



J='(R'^ X M) of closed subsets of M'^ x M is compact [17|, bj and that is 
closed in T{R'^ x 1). 
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Proposition 1. ^usc is compact in the topology induced by lAd- 

Definition 1. Let be (17, ^,P) a complete probability space and let (H, Z) : 

0, — > ^usc be a mapping with 

{uen: T{E,Z)nB^ 0} eA 

for every compact set i? in R'^ x M. Then (H, Z) is called a random marked 
closed set. 

The distribution law of a random closed set is characterized by the prob- 
abilities of hitting compact sets (l8|, ll9(] whereas, by jlSj, Prop. 2.3.1], it 
suffices to restrict to a suitable base. When choosing the same base of all 
finite unions of half cylinders Bi x [tj, oo] as in (23, Thm. XII-6] for random 
u. s.c. functions on W^, we obtain the following characterization of random 
marked closed sets. 

Theorem 1. The distribution of a random marked closed set (as a 

probability measure on ^usc) is completely determined by the joint probabili- 
ties 

P( sup z{x) <ti,BinE^0,ier,BjnE = 0,je{i,...,n}\i), 

where Bi, . . . ,Bn are compact subsets of'R'^, ti, . . . ,tn £ and I is a subset 
of{l,...,n},neN. 

Definition 2. A random marked closed set (H, Z) is called stationary if 

F{T{E,z) + ix,o)e-) = w{T{E,z)e-) 

for all a; G M'^, and it is called isotropic if 

F{eT{E,Z)e-)=¥{T{E,Z)e-) 

for all rotations 9 £ SOa+i with e{R'^ x {0}) = M'* x {0}. 

Example 1. A particular model of a random marked closed set that de- 
scribes an unbiased sampling of a random field [27|| is given when Z is a 
random u. s. c. function on that is independent of the random closed set 
E. We call (S, Z) a random-field model. 

If the data are consistent with a random-field model, any analysis sim- 
plifies considerably since the domain and the marks can be investigated sep- 
arately (see also Remark |4| by using standard techniques for random sets 
[2^ and for geostatistical data 0, For the particular case of marked 
point processes, several tests for the random-field model hypothesis have 
been developed (ill. l23l|. 
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Example 2. Let S be a random closed set and Z{x) = d{dE,x) the Euk- 
lidean distance of rc G M'^ to the boundary of H. Then Z is even continuous 
on H. Since local maxima of Z are only attained at locations in the interior 
of H, the random marked set (H, Z) is a random-field model if and only if 
H = dE almost surely, in which case Z is trivial. 

Example 3. Cressie et al. [3] consider the spatial prediction on a river 
network. Here, H is the flow of the river (as a one-dimensional line or a 
two-dimensional stripe) and Z models the dissolved oxygen. 

Example 4. Let H be a random closed set represented as a locally finite 
union of closed C^-smooth hypersurfaces in such that any two hypersur- 
faces intersect at most in a set of measure zero with respect to the {d — 1)- 
dimensional Hausdorff measure. For any x G E, the mark Z{x) is the maxi- 
mum of the mean curvatures of the hypersurfaces at x. The mean curvature 
has its importance for example in the analysis of foams [l5| . 



3 Characteristics for random marked closed sets 

For the description of random fields a set of second-order characteristics 
like the variogram, the covariance function and the correlation function are 
used 01 . In analogy to these summary functions, several second-order char- 
acteristics for marked point processes have been introduced as conditional 
quantities given the existence of points of the respective unmarked point pro- 



cess [22, l26|]. Since point processes can be described as random (counting) 
measures, these quantities have been derived as Radon-Nikodym derivatives 
of certain second-order moment measures jl, Section 2.7]. Nevertheless, ran- 
dom measures are not always appropriate for the definition of second-order 
characteristics as the following example illustrates. 

Example 5. Let the stationary random closed set H in be given by 
E = C + \J[2z - p,2z + p]U {2z + 1}, 

where p G (0> C uniformly distributed in [0, 1]. Obviously, interpoint 

distances r G (0, 2p] are only possible if both points belong to the same 
segment + [2z — p,2z + p], and interpoint distances r £ {1 — p,l + p] 
are only possible if one point belongs to a segment + [2z — p, 2z + p] and 
the other is from one of the singletons, + 2z — 1} or {.^ -|- 22; -|- 1}. Since 
P(o, r G S) = for all r G {1—p, the approach of defining second-order 
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characteristics using a random measure, which is here based on the Lebesgue 
measure on E}, cannot account for segment-singleton point pairs, and hence, 
these characteristics are undefined for r G (1 — p, 1 Nonetheless, it does 
make sense also to consider the correlation of two marks given that the 
corresponding points are a distance r, r G (1 — p, 1 +p], apart. 

In what follows, B^{x) denotes the Euclidean ball in M"^ with centre 
X ^ and radius e > 0, denotes Minkowski addition, and we write 
shortly for H © Be{o). 

Let (H, Z) be a stationary random marked closed set in with marks 
in M. For any £ > define the (stationary) random field by 



Let / : ^ M be a right-continuous function. For all h G define 



whenever i^\f\{h) < oo and P(o, /i G "©e) > for all £ > 0, otherwise Kf{h) 
is undefined. 

In particular, for the following choices of /, 




otherwise. 




(2) 



e(mi, m2) = mi, ci 



'{mi,m.2) = mim2, v{mi,m2)=mi, 



(3) 



define 



E{h) = Ke{h) 

7(^) = ^(/^f (M + Ky{-h)) - Kc{h) 

COv{h) = Kc{h) — Ke{h)Ke{ — h) 



(4) 



(5) 
(6) 



Kc{h) - Ke{h)K,e{-h) 



(7) 



kmm{h) = (m)"^ Kcih), (m 7^ 0), 



(8) 



where 



m = E[Z{o) I o G S] 



is the mean mark. 
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We call 7 the mark variogram, cov the mark covariance function, cor the 
mark correlation function and kmm Stoyan's kmm-function of |22]. 
Note that, if H = R"^, these definitions are compatible with the classical 
definitions for random fields (see Remark |4]). 

Whenever (S, Z) is assumed to be both stationary and isotropic the char- 
acteristics given by (I4l)-(l8|) are rotation invariant. By slight abuse of notation 
we will write E{r), r G [0,oo), instead of E{h), h G M*^. The same applies 
for the functions defined in Eq. ©-Jll). 

Remark 1. Let = -nH^^) be the random volume measure associated 
with the random closed set Here, is the d-dimensional Lebesgue 

(2) 

measure. If //g denotes the second-order moment measure of ^'g then, for 



dx dy 



Bi,B2e BiR"^), we have 



E 
E , , 

' 'B2 JBj 



B1XB2 



B2 JBi 



f[Z,{x),ZM) ^e{dx)^>e{dy) 
/(mi, 7712) Qe;x,y{d{mi,m2)) ^if\d{x, y)) 
fimi,m2) Qe-x,y{d{mi,m2))Fix,y G H®£)dxdy, 



where Qe;x,y is the two-point mark distribution of the weighted random mea- 
sure (^'g,^^) Hence, for almost all {x,y) with P{x,y G S^e) > 0, we 
have 



E 



/ \Ze{x),Ze{y)] I x,y G 



/(mi,m2) Qe-x,y{d{mi,m2)). 



Remark 2. In case P(o, /i G H) > 0, /i G W^, the above definition takes the 
simpler form 

Kf{h)=nf{Z{o),Z{h)) I 0,/iG H] 
if we impose the integrability conditions 

E[|Z(o)|1h(o)1h(/i)] < 00, K|g|(/i) < 00, 

for / = e, and, 

E[|Z(o)pls(o)lH(/i)] < 00, K\^\{h) < 00, 
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for / = c or / = v. This can be seen as follows. Denoting by a+ and a_ the 
positive and the negative part of a G M, respectively, we always have 

^e(o)+lsee(o)lsee W < ^^(o)+lse.(o)lse.(^) < \Z^{o)\l^^,{o)l^^,{h) 

for all e < e, where the right-hand side is integrable due to n^f,^{h) < oo. 
Similarly, 

Z,{o).ls^^io)ls^M < Z{o).Uio)lsih) < \Zio)\ls{o)lsih). 
In the same way we obtain 

\Z^{o)\H^^,{o)l.^,{h) + \Z{o)\'U{o)l.{h) 

as an integrable upper bound of |Z£(o)pls^^(o)ls^^(^)| and 

{\Z{o)\' + \Z{h)\') ls{o)ls{h) + (|Z^(o)|2 + \Z^{h)\') lse,(o)lse. W 

as an integrable upper bound of \Z^{o)Z^{h)\l'^^^{o)lw^^{h). Since Z is 
u. s. c. on H, a value e > exists for every x G E and for every 5 > such 
that Z{y) < Z{x) + S for all y G B^ix) n E. Hence, we have ^^(x) Z{x) 
from above as e — > 0+. Further, x ^ S implies x ^ for all sufficiently 
small e. We then have 

f{Uo)X{h))UeM^^^M ^ f{Z{o),Z{h))U{o)U{^ a.s. 

as e — > 0+. Hence, by the dominated convergence theorem, we have 



i^f{h) = lim 



E 



/ (z,(o),Z,(/i)) lsee(o)lseeW 



£-0+ P(o, h G See) 

^ E[/(^(o),Z(/^))ls(o)l5W] 

P(o, /i e s) 

Remark 3. There exists an alternative concept of random marked sets which 
is inspired by the notion of random fields and where second-order character- 
istics in the sense of the preceding remark can be defined. 

Let M0 = R U {C0} be the extension of R by some C,0. We denote by 
i3(R0) the respective Borel cr-field which is generated by all sets Bi U B2 for 
Bi G ^(R) and B2 C {-00,00,(0}. 

A family of random variables Z{-,x) : — > R0, x G R'^, on the probabil- 
ity space (J7,.4, P) is called a random field with random domain S, if 

S = {xGR^: Z(-,x) /C^}- 
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Clearly, when Z takes only values different from C0 or {—00,00,(^0} this 
notion of a random marked set includes usual M- or M- valued random fields 
on W^. 

Note that H is a random set in a very general sense entirely deter- 
mined by its indicator Isix) = l^{Z{x)). If Z is jointly measurable, i.e., 
Z is {A0 i3(M'^), i3(M0))-measurable, then the realizations of S are almost 
surely Borel measurable. If we have even almost surely closed (open) realiza- 
tions of H then Z is called a random field with random dosed (open) domain, 
see also [l9| . 

If P(o G H) > holds for a stationary random field Z with random 
domain H we can define second-order characteristics without any further 
assumption on path regularity. Let Z be the (stationary) random field given 
by Z{x) = Z{x) for x G S, and Z{x) = otherwise. Let / : ^ M be a 
measurable function. For all h define 

Kf{h)=¥.[f{Z{o),Z{h)) I 0,/iG H] 

whenever P(o, /i E S) > and E[|/(Z(o), Z {h))\l^{o)l^{h)] < 00. 
Remark 4. Let (S, Z) be a stationary real- valued random-field model and 

max Z{y), x £ H^e, 
ZAx) = i yeHnB,(x) 

Z{x), otherwise. 

Since Z is u.s.c. on H we have Z^{x) Z(x) from above for x € H, and 
hence, by the definition of Z^, for all x G M*^ as e ^ 0-|-. Then, using the 
independence of Z and H, we obtain 

AJ/(/i) =^hm^E[/(Z,(o),Z,(/i)) I 

= hm E[/(Z,(o),Z,(/i))] 
e^O-l- 

= E[/(Z(o),Z(/i))] 

for all /i G R*^ which satisfy P(o, h G H^^) > for all e > and, depending 
on the choice of / according to ([3|), one of the integrability conditions in 
Remark [2] with H replaced by . 

Remark 5. The definition of according to ([2]) is, in important situtations, 
consistent with the classical definition of the second-order characteristics of 
stationary marked point processes [23|: Let $ be a stationary simple marked 
point process on M'^ x M. Then, H is the support of the unmarked point 
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process <I> = <!>(• x M). We assume that the second-order moment measure 
^(2) Qf (J) jg locally finite. For \\h\\ > we have 



E 



/ [Ze{o),Zs{h) j lsS)B4o){o)'\-S(BB,{o){h) 
= E / (^Ze(o), li^^(^Be{o))=l}'^{<S>{B4h))=l} 
+ E f (^Zs{o),Zs{h)^ l{<I>(B4o))>l}l{<J.{B^{/i))>l} 
+ E f iZs{o),Zsih)] l{<I>(B,{o))=l}l{${B^{h))>l} 



For any < e < ||/i||/2 the first summand equals 



E 



^ f{mi,m2)lB,{o){xi)'i-B,{h){x2) 



We can extend the argumentation in [8|, Prop. 9. 3. XV] in order to conclude 
that 



mB,io))>l,^B,ih))>l) 



1 



/i(2)(S^(o) X B,{h)) /i(2)(S^(o) X 

as £ ^ 0+. If we additionally impose the condition that for some e > 0, 



E 


/(z,(o),Ze(/l)) 


l{*(i3ao))>l}l{*(S.(/i))>l}l{|/(z,(o),Z,(h))|>A/} 


P 


mB,{o))>l,^B,{h))>l) 



sup 

e6(0,e) 

as M ^ oo, we obtain 



flf\B,{o) X Be{h)) 



Kf{h) = lim 



0+ ^^(^){B,{o)xB,{h)) 



which equals ^^^^-a. e. the Radon-Nikodym derivative 

d^j (j;, + h) 
d^(2)(a;, a; + /i) 

For instance, the above condition is satisfied if E 
formly bounded on (0,e) for some a > 1. 



f(Z,{o),Z,{h) 



IS uni- 



10 



A function / : R'^ ^ M is called positive definite if 

n n 

aiajf{xi - Xj) > 

for any n G N, . . . , x„ G M*^, and ai, . . . , a„ EM, and / is called condi- 
tionally negative definite if 



aiajf{xi - Xj) < 

i=i i=i 

for any n G N, xi, . . . , x„ G M'^, and for all oi, . . . , a„ G M with flj = 0. 

For a random-field model all second-order characteristics coincide with 
that of a random field with u. s. c. paths, see Remark HJ and thus, share the 
same definiteness properties. On the other hand, for marked point processes 
it has been shown by examples 28 1 and systematically that the mark co- 
variance function, the mark correlation function and the A;mm-function need 
not be positive definite, and the mark variogram need not be conditionally 
ne gati ve definite in contrast to random fields. Some of the constructions used 
in [22] are based on the fact that for a marked point process, H is a locally 
finite subset of and has therefore Lebesgue measure zero. However, the 
next example shows that in general we cannot expect that the mark covari- 
ance function is positive definite (and the mark correlation function and the 
^mm-function either) even when we have E[frf(H n [0, 1]*^)] = P{o G H) > 0. 



Example 6 (Continuation of Example [2]). Let be p G (§,!], ^ a random 
variable uniformly distributed on [0,1], H = 1 
1-periodic function defined by 



Zi^,x) 



^ © [C,P + C], and Z{^,-) a 

P-{x-0, xGZe[£ + e,p + 0, 
0, X G ze + + 



Then Z and H are jointly stationary and each of the characteristics given by 
©-([HI) is 1-periodic. In particular, on [0,1/2) we have 



r G [0,1-p), 



GOV r 



' p4_4p3y_;^2p^r^+48pr3-36r'^ 
48{p-r)2 ' 
-3p"'-24p'^7-+12p^-48p^r^+48p^r-12p^ 
48(2p-l)2 

I 64pr3+48pr2-48pr+8p-32r3-24r2+24r-4 ^ r-, p\ 

48{2p-l)^ ' ^^i^-P^2> 

4p"'-8p^+6p^-2p+12r-^-24r3+18r^-6r+l r a \P i] 
12(2p-l)2 ' ^'=i2'2J' 
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and, by symmetry, cov(r) = cov(l — r) for r G (^, 1). Since cov is 1-periodic 
the 0th coefficient of the Fourier series of cov is proportional to 



1 

J cov(r) dr = In 



_l_ 409p'^-790p^+565p3-280p^+120p-24 



^2p-l J ' 120{2p-l) 



which is negative for | < p < 1 (and vanishes for p = 1, which is the random 



field case). Since cov is continuous, Bochner's theorem [2li| implies that 
cov(r) cannot be a positive definite function. 

Our major example is analysed within an own section since some results 
might be of interest not only to the field of random marked sets but also to 
the theory of positive definite functions. 



4 Gaussian random fields exceeding t eM. 

Let Z he a stationary and isotropic centered unit variance Gaussian random 
field in W^. Then, for t e M, we define 

Et = {xeR'^ : Z{x) > t}. 

If, in particular, Z is almost surely continuous [H, Q] then is almost surely 
closed, i. e., (Ht, Z) is a random marked closed set. Note that Ht is a so-called 
excursion set which has been extensively studied in the literature, see [H, [3] 
and the references therein. 

Since Z is assumed to be both stationary and isotropic, its covariance 
function, Cov(j;,y) = E[Z(x)Z(y)], x,y £ M."^, is translation and rotation 
invariant, i. e., there exists a function R : [0, 00) R such that 

Cov{x,y) = R{\\x-y\\), x,y eM.'^. 

First, we consider the case t = 0. 

Theorem 2. Let Z be a stationary and isotropic centered unit variance 
Gaussian random field in M*^ with covariance function given by R : [0, 00) 
M. Then, forr G [0, 00), the second-order characteristics of{Eo,Z) are given 
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by 



E{r) = ^±^^^ (9) 
V 2 arcsin(ii(r)) + f ' 

, . N Jl-RirY vr {l+R(r)f , , 



, . ^ fi(r)(arcsin(fl(r))+|)2 + 0^^1(O^(arcsin(fl(r))+|)-f(l+fl(r))^ 

^ (arcsm(i?(r))+§)2+_R(r)7l-iJ(r)2(arcsm(i?(r))+§)-f (l+_R(r))2 ' ' 

Obviously, each of the second-order characteristics of (Hq,^) is a con- 
tinuous transform of R. In particular, this means that continuity of R is 
preserved. Vice versa, due to the monotonicity of the transform for cov (see 
Theorem [3] below) already from cov it can be deduced whether or not R 
contains a nugget effect [6||. 

Since, for every stationary random-field model, E{h), h G M'^, is constant, 
equality ([9]) implies that (Hq, Z) is not a random-field model, i. e., there does 
not exist a random field in whose second-order characteristics coincide 
with that of (Hq, Z). It is therefore quite surprising to see that we are not able 
to falsify that (Hq, Z) is a random-field model by using the mark covariance 
function or the mark correlation function of (Hq, Z). 

Theorem 3. The functions /o : [—1, 1] M, 



Up) = p + 



arcsinp-l-^ ^ (arcsinp -|- 
and go : [—1, 1] M, 



^ ^ p(arcsinp-F ^)'^ + -Jl - p"^ {avcsin p + f ) - ^(1 + p)'^ 

9o{P) = ^ , 

(arcsin/9 -|- |)2 -|- pyl — ^^(arcsin p + f ) — f (1 + p)^ 

are absolutely monotone on [0, 1], i. e., they have only nonnegative derivatives 
there. 

Corollary 1. cov(r) and cor(r) are positive definite functions. 
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However, the following proposition shows that the /cmm-function of (Hq, Z) 
is not positive definite, in general. 

Proposition 2. Let Z be a stationary and isotropic centered unit variance 
Gaussian random field in with covariance function given by a continuous 
function R : [0,oo) — > M. Then the kmm-function of {r.Q,Z) is positive 
definite if and only if R= 1 . 

Proof. Let 



QiP) = P + 



arcsin(p) + ^ 
Then q{l) = 1 and 



q{p)>p + 7:-Wl-p^ > 1, PG ((7r2-l)/(7r2 + l),l). 

Hence, kmm{f) is not a positive definite function if i? ^ 1 Theorem 
1.4.1]. □ 

The mark variogram of (Hq, Z) is in general not conditionally negative 
definite, which can be seen as follows. Consider (Hq, Z) for dimension d = 1 
and R{r) = cos(r). Since 7 is conditionally negative definite if and only if 



e"'''^ is positive definite for all s > 2l|, Theorem 6.1.9] it suffices to show 
that e"*^ is not positive definite, e"'''^''^ inherits 27r-periodicity from cos(r), 
and hence, it is positive definite if and only if its Fourier coefficients are 
nonnegative. Numerical calculations yield that the first Fourier coefficient is 
nearby -0.03364. 

Now we switch over to the more general case t G M. Unfortunately, 
unlike the case f = 0, we cannot express all the second-order characteristics 
of (Hf, Z) in closed form. In particular, for a stationary and isotropic centered 
unit variance Gaussian random field Z in with covariance function given 
by R : [0, 00) M, we have 



RUM) 



F{o,heEt) 
see [i, Eqn. (10.8.3)]. Here, 



j ip{t,t,s)ds + ^{tf, 



^(*) = / ds, <^{t) = —=e 2 , te 



t 
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denotes the tail probability function of the standard Gaussian distribution. 
By ip{s, t, p) we denote the density of the bivariate Gaussian distribution 
with unit variance and correlation p. In the following we concentrate on the 
mark covariance function of {Et,Z). 

Lemma 1. Let Z be a stationary and isotropic centered unit variance Gaus- 
sian random field in with covariance function given by R : [0, oo) — > M. 
Then, for i G R and h G R'^, we have 



E[Z(o)l[t,o,)(o)l[t,„o)(/i)] = Etm\h\\)), Et{p) = ip{t){p + l)^i^t^j 
and 

E[Z(o)Z(/j)l[,,^)(o)l[,,^)(/i)] = CtiRiWhW)) 



l-p 



+ P 



where 



Ct{p) = r ip{t,t,s){p-s-t^)ds+p^\t)+ip^{t) 

Jo 



If we write 



Ptip)= [%it,t,s)ds + ^{tf, pG[-l,l], 
Jo 



then the mark covariance function of {Et,Z) is given by 

cov(r) = ft{R{r)), 

where 

Ctip) Mpf 



hip) 



Ptip) Pt{pf 



There is strong evidence that also the mark covariance function of {Et,Z), 
t 7^ 0, is positive definite for a certain class of Gaussian random fields Z. 
Figure m shows ft{p) and fl{p) for several t, indicating that for these t the 
functions ft{p) are both increasing and convex for p G [0,1]. Hence, if this 
is really true, for instance Polya's criterion would imply that, for any 
continuous and convex function R : [0, oo) M satisfying R{0) = 1 and 
lim^^oo = 0, the function ft{R{\ • |)) is positive definite on M. 
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Figure 1: ft (left) and // (right) for i = — 1 (dashed), t = (dotted) and 
t = 1 (solid). 



5 Differentiability at 

In this section we continue the example of the preceding section and show 
that the corresponding mark covariance function has a right-hand derivative 
at which does not vanish. 

Lemma 2. Let Z be a stationary and isotropic centered unit variance Gaus- 
sian random field in with continuous covariance function given by R : 
[0,00) excluding R=l. Let be Chj(||/i||) = ¥{o,h £ Ef) the set covari- 

ance of the excursion set H^. If Z is mean-square differentiable then 

ci(o+) = -^V-^"(o+)<o, 

otherwise (74^(0+) = —00. 

Theorem 4. Let Z be a stationary and isotropic centered unit variance 
Gaussian random field in with continuous covariance function given by 
R : [0,00) M. excluding R = 1. Then the mark covariance function 
cov(r) = ft{R{r)) of (Ht,Z) has a negative right-hand derivative at r = 0, 
in particular, we have 



{t^-l)^itf-3t^m{t) + 2^{t)' ^(t)y/-R"{0+) 
in case Z is mean-square differentiable, and cov'(O-l-) = —00 otherwise. 
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6 Proofs 

First, we prove Lemma[T]of Section [4] which is needed in the proof of Theorem 

El 



6.1 Proof of Lemma [T] 

The proof is based on techniques used in [3, p. 351]. We have 



oo poo 

Et{p) = I J xnp{xi,X2,p)dxidx2 

foo ;>oo / 1 />oo POD 



Xl 



t Jt 

fOO fOO 



(47r2 J_ 



oo J — oo 



4^2 / / Eye-it(yi+s/2)_ I _ _ _ ) dys 



CXD J — OO 



1 / t 1 

W2 Viyi ^1 



after reversing the order of integration and integrating out xi and X2. For 
the part 



I rcc. rco it{yi+y2) 

EtAp) = ^T^ / e-^^^ -^dyid2/2 

J —OO J —CO 



-viy2 



we have 



JO 



For the part 



Et;2{p) 



1 

4^ 



oo roo 



e 2 



-it{yi+y2) 



oo — oo 



dyi dy2 



we obtain 



d''Et.,2{p) i 



47r2 



oo roo 



y2e 



— oo J — oo 



dip{xi,X2,p) 



dx-} 



(f{t,t,p) 



dyi dy2 
t 



{xi,X2) = it,t) 



By Et{0) = (p{t)'i/{t) and 



dEtip) 



dp 



/ Xl - 



dp 



dxidx2 = {tip{t) + ^{t))ip{t), 



p=0 
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we get 



Et{p) = t I I ^l;"';';^ dsdp + t £ if it, t, s) ds + {p + 



D Jo 
(P+1) 



Jo 1 + s 



Jo 1 



+ s 



ds + (^(i)*(i) 



Further, 

p 



1 



l + s 



-fu 



27r 7 ^/2u-l 
1 

tV2 



d-u = 



7 



47r 7 -y/u 



(14) 



27r 



/ 



ty/{i-p)/{i+p) 



which gives the first equahty in the lemma. 
The second equahty is shown similarily: 

Ct{p) = E[Z{o)Z{r)^,oo){o)l[t,oc){r)] 

XiX2(p{xi,X2, p) dxi dX2 



t Jt 

oo fOO 



XiX2 



(47r2 



oo /•oo 



oo J — oo 



T 













dj/i dy2 dxi da;2 



1 



-j^ /"OO /"O" 

^ ^ / 

J — oo •/ — oo 

Thus, we can write Ct{p) = Ct-\{p) + Ct-2{p) + Ct-3{p), where 

r e-i*(2/l+y2) 

Q-^y dyidy2 



dyi dy2. 



-| /" OO /" cx; 



oo ./ — oo 



-yiy2 



CtM=tj^2 



oo />oo 



e 2 



-it{j/i+y2) 



1 



oo J — oo 
oo /•oo 



47r2 



oo J — oo 



-it{j/i+'y2) 

— 2~ ^2/1 dy2 

-12/1^2 
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and 



' — OO J —CO 

For the latter expression we have 

47r2 
= Ptip) 



— "iyi ^2^2 



5p 



e 2 



-it(yi+y2) 



OO ^ — OO 



Since Ct(0) = (p(t)l\J2Tx we then have 

G(p)=2t2 r rf^MA^sdp+ r f\{t,t,s)dsdp 

Jo Jo ^ + s Jo Jo 

fp 

+ / ip{t, t, s) ds + p^{tf + 2ptip{t)'^{t) + i/{t). 
Jo 

fp 

= / ^{t,t,s){p- s-t^)ds + p^'^{t) + ^'^{t) 
Jo 

+ 2t^{t) \^{t) + [p+ 1)^ (t^{i-p)/(i + p) 



6.2 Proof of Theorem d 

As 

P(Z(o) > 0, Z{r) > 0) = ^ (arcsin(i?(r)) + |) 
the formulae for E and cov follow immediately from Lemma [H Now, 

vr ' 

and it remains to show that 

E[(Z(o))2lp,^)(Z(o))l[o,oo)(^(r))] 

_ R{r)^Jl- R{rY + arcsin(i?(r)) + f 



m = E[Z(o) I Z(o) > 0] 



27r 



We have 
Voip) 



OO /•OO 



x\ (f{xi,X2,p) dxi dX2 



JO 

OO roo 





2 
4^ 



(4vr2 



OO />oo 



xje-'y "^e-^y ^ydyidy2 ] dxidx2 



OO J — OO 



OO /"OO 



OO J — OO 



"1— dyi d7/2- 
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The coefficient of p in the exponent of the integrand is (—2/12/2)- Hence, 
differentiating the last equation three times with respect to p yields 

^^^LlJ' ' '''''' 

27r 27rVdet(S-i) loo i-oo ' 

_ -2Vdet(S-i) 1^ _-2 1 



27r ^ ^22 27r _ ^2 ^ 



Because of 



dp^ 47r2 yi 
we have %^(0) = 0. Hence, we obtain 

9Vo(p) 



,2 



after integration with respect to p. Because of ^-(0) = l/vr and Vo(0) = 1/4 
we then obtain after further integrations with respect to p 



yo{p) = ^ [P\/ 1 - p2 + arcsin p + - 



6.3 Proof of Theorem H 

The idea of the proof is to show that the coefficients of the Taylor expansion 
of /o and qq are all nonnegative. We give the proof for /o only; for qq, the 
same techniques can be applied. 

Since arcsin can be extended to the complex plain C by 



arcsin z = — iln(i2: + \/l — z^)^ z G C, 

with branch cuts (— 00,— 1) and (l,oo), also /o can be extended to C. We 
restrict our attention to 



arcsin z + ^ ^ (arcsin z + 

Since arcsin z = — 7r/2 is satisfied only for z = — 1 and the denumerators of 
all derivatives of / consist of a product of integral powers of arcsin z + 7r/2 
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and Vl — z^, / is not only holomorphic in the open unit disc D = {z E 
C : |2;| < 1} but also on C \ ((— oo, — 1] U [1, cxd)). This holds because -y/i is 
holomorphic everywhere except on ( — oo, 0] where ^/z is not even continuous. 
The expansion of arcsin z around z = ~1 

arcsinz = -f + + z{l + 0(1 + z)) , \z\ < 1, 

shows that the derivatives of any order of / can be continuously extended 
in 2; = — 1 as z — > —1, \z\ < 1. Hence, the only branching point of / in D is 
z = 1. From the expansion of / around z = 1 

/W = -f +Cl(l-z)U(f + (l-^)+C2(l-^)i 

+ 0{{l-zf), \z\<l 

with 

ci = -V2(^-i), c, = V2(^ + ^ + ^-i|), 

we see that the asymptotic behaviour of / around z = 1 is essentially de- 
scribed by 

g{z) = ci{l - z)^ + C2{1 - z)^ . 

At first, we derive an upper bound for the coefficients bn of the Taylor ex- 
pansion of 

h{z) = f{z)-g{z) 

at 2; = 0. The function h is holomorphic in D and by construction twice 

continuously differentiable in D since h^'^\z) behaves around z = 1 like 
1 

c + 0((1 — z)2), |2;| < 1. As a consequence of Cauchy's integral formula we 
have 

for all £ G (0,1) and -Di-e = {2; G C : < 1 — e}. After twice partially 
integrating we have 

1 1 f h^^HC) 

K = 7r^-l TT / T^^^' ^ = 2,3,4,.... 

27ri n{n — 1) J ^ 

dDi-e 
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Hence, for n > 2, we obtain 



l^nl < r(l - e)^"" max 

2tt n{n — 1) 'p&[o,2tt] 



/i(2)((l-e)e'^) .27r(l-e). 



is continuous in D we can conclude that for n > 2 



\bn\ < 



max 



n(n — 1) ¥'G[0,27r] 



/i(2)(ei^) 



< 



1 



n{n — 1) 



0.182 



by Lemma [3] below. 

In the following we derive a lower bound for the coefficients o„ of the 
Taylor expansion of g{z) at z = 0. Note that we have the series expansions 



1 °° /1\ °° 

n=0 ^ ^ n=l 



(2n - 2)! 



n!(n - l)!22'i 



-z", \z\ < 1, 



and 



(i-^)t=i-i.+3f: ..,p":l_, .". N<L 



n=2 



n!(n-2)!22"-2' 



Hence, for n > 2, we have 



(2n-2)! (2n-4)! 
«n = -Cl • — + C2 • 3 ■ 



n!(n- 1)!22 
Since ci < and C2 > 0, we obtain 
(2n - 2)! 



(n-2)!22"-2- 



On > -Cl 



n!(n - l)!22"-i - ^^2/^ y ^ 
by Stirling's approximation. At least for n > 6 we have 



e / n — 1 \ „3 



2e 12n 12(n-l) 



1 + ^ ) ^ 12{n^l) > M + ^ 1 e > 1. 



1 



Thus, 



a„ > 



4 - vr 
7r2-v/27r 



n 2 . 



n > 6. 



Combining this with the upper bound for it is easily seen that the nth 
coefficient of the Taylor expansion of /o at x = is nonnegative for all n > 30. 
The remaining coefficients of order n = 0, 1, . . . , 29 can be determined from 
a series expansion and are easily evaluated to be nonnegative. 
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Lemma 3. We have 



max 

(pe[0,27r] 



< 0.182. 



Proof. Figure [2] indicates that this maximum is achieved for f = 0. From 
the expansion of / around z = 1 we obtain the value of this maximum: It is 
twice the absolute value of the coefficent of (1 — z)^, i. e., we have 



max 

VG[0,27r] 



h('\e^n = \h^'Hl)\ = {^-^ + ^ + ^-^)< 0.08. 



Formally, we can actually prove that an upper bound is 0.182. Since the 
conjugate- complex of /i^^^(e"^) equals /i^^^(e~"^) we can restrict to the interval 
[0,7r]. 

First we make an estimate for [f,vr] then for [0,^]. From the series 
expansion 

n=l ^ ' 

we obtain the two series expansions at z = — 1 

1 ^ / ^T^_./7t/^, (2n-2)! (\^z 

n!(n- l)!22""i 



:i-.)2 = v2^i-i2Li = v2(i-5; 

oo 



n=0 

and 



arcsinz + f = -IML^ ^^(i + ^)-+^ 

2 ^ n!n!23"+i n + ^ 



oo 



(l + z)5^n«(l+zr, 



n=0 

which are convergent for \\^z\ < 2. From the latter we also obtain according 
series expansions 

oo 

(arcsin z + f = (1 + t ^ u^f ) (1 + z)", A; = 2, 3, 4. 

n=0 
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Let 

(z) = (1 + z) 1 J] ) (1 + z)", A; = 1, 2, 3, 4 



N 



n=0 



(u\ , , arcsin z + ^ — Pil?^ iz) 
V%\z) = , \r ^ % ^ = 1,2,3,4, 

{1+ Z)2 



and 



(l-z)2 -Q^(z) 



n=0 



(1 + ^) 



We write /i*-^-* (z) as a fraction with numerator 

hnum{z) = 2(arcsinz+ + z){l - z) - (arcsin 2; + f )^(1 + z)(l - z) 



- 7r(arcsin 2; + f )^(1 + 2) 5 (1 - S 



+ 2;(arcsin2; + |) (1 + 2)2 (1 - 2)2 - z"'(arcsinz + f ) 
+ ^(arcsinz + f )^(1 - 37r(l +z)t(l - z)^ 



+ 47r (arcsin 2; + f)(l + z)^(l - z) 

+ 7rz(arcsinz+ f)(l + - ^(arcsinz + f )^(1 + 2;)i(l - 

and denumerator 

^denum(^) = (arcsin z + f )*(! + § (1 - z)i. 

By /tnum°^(2;) we denote the approximation of hrmmiz) if we replace (arcsin z+ 
f by Pi}\z) and ^/T~l by QaK^) within /inuni(-z)- On the set {e^*' : (p G 
[^,7r]} we can make the following estimate for the error 



Since we have 



^denum (s"^) 



|(arcsinz + f)^| > 4|l + z| 
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we obtain 

r,{cp) < 1|1 - e^^|-i (2|pSi)(e'^)| • |1 + e^^| • |1 - e^'^| 

+ |pg)(e^'^)|-|l + e''^|-|l-e''^| 

+ 7r|pg)(e^'^)| . |l + ei'^|i • |1 -e^'^l • |Qiv(e^'^)| 
+ TT\pg\e'^)\\l+e'^\l\l-e'^\-\qN{e'^)\ 

+ \p(-^\e'^)\.\l + e''^\"^-\QNie'n\ 

+ iPj^He'ni -11 + e'l^ • kiv(e^'^)| + bg^(e'^)| • |1 + e'^l^ • |gjv(e^'^)| 
+ IP^i^'ni + ^l#(e''^)l • |1 + e''^|i +37r|l + e''^|i ■ \qM{e'n\ 
+ 47r|l + e'^|2 . |1 - e''^! • \p^^\e''^)\ + 7r|l + e^'^l^ • \p^^\e''^)\ 

+ ^|pS)(e-)|.|l + e-|-|l-e-|). 

Since |1 + e"''| is decreasing on [0, tt] we have 

ll + e^'^l < |l + e^i| = VS 
(k) 

on [f ,7r]. The coefficients Un , k = 1,2,3,4, and Vn are non-alternating for 
n > 1. Hence, on [§,7r], we have 

\p^^\e'n\<\p'^^\e'n=:d'^^\ k = 1,2,3,4, 
kjv(e^^)|<kiv(e^^)|=:eiv 

for all A?^ > 1. Furthermore, on [f ,7r], we have 

|l_e'^|-i < 1 
and, for A'' > 20, we can estimate 

|Qg)(e'^)| < 1.001- |l-e'^|i 
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Hence we obtain 



ri < i(2V3(5{^^ + VsdP +^33^^^ • 1.001 + 2.12^3 



3 

4 EAT 



+ vr 6^ EN + 3i<5g^ • 1.001 + 2.122 . 3!^^)^^ + §^ 



(3) 

N 

3 -(4)^ 



^ ■ ' iv ' iv ' ^ 

Thus, for > 20 we have rj < 0.005041. Furthermore, we can estimate 



T approx 
"-num 




< 


.approx / i<^^ 
"■num ) 


^dcnum 




7 

4|1 + e^f\2 



By construction, the right-hand side of this inequahty is the absolute value 
of a polynomial in 1 + e"^ which takes its maximum on , vr] at = tt. For 
= 20 this maximum value is less than 0.1764. Combining both estimates 
we end up with max^^^yr |/i(2)(e"^)| < 0.182. 



Similarly, from ([15]) we obtain at z = 1 the series expansions 

(2n-2)! 



00 



l-z 



^ n!(n - l)!22"-i V 2 



1 - z 



n=0 



{l + z)-2=V2j2 



(2n)! 



n=0 



n!n!23'^+i 



(i-zY 



and 



n=0 



arcsin z + 2 = tt 



^E 



(2n)! 



n=0 



which are convergent for \l — z\ < 2. From the latter we also obtain according 
series expansions 

00 

(arcsin z + ^)'' = tt'' + ^ u^J^\l - z)f , k = 2,3,4. 

m=l 
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Let 

Pff)(z)=7r^ + |;SL^)(1-Z)f, fc = 1,2,3,4, 

m=l 

N 

n=0 

N 

RN{z) = J2wn{l-zr 



2N 



n=0 



and 



ti.) , , arcsin z + 5 — plf -* iz) 

Vn\^) = 3 ^ % ^ = 1,2,3,4, 

(1-z)- 



„12 

1 



(In{z) 
rN{z) 



- z 



{l + z)2-Q^{z) 

(i-z)i ' 

(1 + - Rn{z) 
(l-z)i 

We write ^^^^ (z) as a fraction with numerator 

hnnm{z) = (l + z) 2 (2(arcsin z + f )(1 + 

- (arcsin z + + z)(l - z) 

- 7r(arcsinz + f )^(1 + 2)^(1 - z)^ 

+ z(arcsinz + + z)^l - z)^ - z^(arcsin2; + 

+ J (arcsin z + f )^(1 + z)t - 37r(l + z)^l - z)^ 
+ 47r(arcsinz + f )(1 + zf{l - z) 

+ 7r2;(arcsin2; + f )(1 + zf - ^ (arcsin z + + z)i(l - z)) 

and denumerator 

^denum(2;) = (arcsin z + 1)^(1 - 

By /inum°^(-z) we denote the approximation of /inum(^) if we replace (arcsin z+ 
If by P^^\z), by Qiv(z) and (1 + z)-\ by i?iv(z) within K^r.{z). 
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On the set {e'*' : <p E [0, 5]} we can make the following estimate for the error 



^denum (g'''') 



Since on [0, |] we have 



I arcsin(e"^) + ^| > | arcsin(e'3 ) + ^| > 2.11 



we obtain 
2.11%ip) < \KumieV\rN{e'n 



+ 



2|#(e'^) 

5(2) 



11 + e 



e'1 + |pS^(e'^)||l + e 



1(^1 



+ 7r|P^^^(e^'^)||gjv(e"^)||l + e^'^||l-e'^|2 
+ Tr\j^^\e'^)\\QN{o'n\\'^ + e'^\\l-e'^fi 



+ 7r|p<^)(e'^)||g^(e'^)||l + e 



+ m|pi;J(ei'^)||Q^(ei'^)||l + ei'^| 



+ ^|p)^^(e'^)ll?iv(e^'^)||l-e'^|2|l + e^^| 

+ 37r|giv(e''^)||l-e''^|^|l + e'<^|' 

+ 47r|l - e'<^||l + e''^f\^^\e''^)\ + 7r|l + e''^|2|p<J(e''^) 

3|C2| 



4 

3|C2||~{4) 



^'F(e^'')ll9iv(e'^)||l + e^'^||l-e" 
l#(e''^)IIQiv(e^'^)||l+e^'^||l-ei'^| 



31^21,^4) 



|p)^'(e^'')ll5iv(e^'^)||l + e'n|l-e^'^|2 
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Since |1 — e"''| is increasing on [0,7r] we have 

ll-e'^r^ < |l-e^i|-^ = 1 

on [0,^]. The coefficients Un \ Vn, and Wn are non-alternating for n > 1. 

Furthermore, for m > 1, we have > *-* l^2ml > l^2m-i-il' ^ ~ 2)3,4. 
Since we have 

|l-e'^| < |l-e^i| = 1 
on [0, |] we can conclude that there 

|p!;^(e^'^)l < |pS^(e^^)l =: 4'^ k = 1,2,3,4, 
\qN{e'n\<\QN{e'^)\=-^N, 
\rN{e'^)\<\rN{e'^=:KN 
is satisfied. Finally, we have 

|Pif^(e"^)| <^', fc = 1,2,3,4, 
|l + e''^|<2, |Q7v(e''^) < 2. 



Hence we obtain 

r? < (47r + 2Tr^ + 2V27r^ + V2tt'^ + tt^' + ^\/27r^ 
+ 12\/27r + levr^ + 47r2 + ^V2Tr^)KN 
+ 255^)3-3 + 5^^3-2 + tt^enS-^ + TT^S^^ + 7r5S^^£Ar3-5 
+ Tr^ejv + ^23"^ + 5g)£^ 

+ + M^4~^3-i + ^ M44)e^3-i + 37r?^ V2 + 4;rv^« 

+ 7rV2^g) + f |c2|7r^?'iv3-i + ||c2|4' + ||c2|4'?'iv3-^ • 
For iV > 12 we have rj < 0.01052. 

Both the functions | arcsin(e''^) + f | and |/inS?n™(e''^)(l-e'^)-3/2l for TV = 12 
are decreasing on [0, In order to estimate the quotient of both we divide 
[0, |] into small subintervals of length For each subinterval it is easy to 
evaluate that the quotient is less than 0.1. Hence we have 

Tapprox / 

^— ^ < 0.1 

^denum (s^''') 
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n 1 



5 6 



Figure 2: \h^^\e''^)\ for ^ E [0,27r]. 



which imphes together with the estimate for r] that also on [0, ^] |/i'^^)(e"^)| 
is less than 0.182. □ 



6.4 Proof of Lemma [2] 

We have -Pt(l) = P(Z(o) > t) = ^'(t), and a series expansion of Pt{p) around 
p = 1 yields 

Pt{p) = ^{t) + ^ (-2yw + o{{i - p)!; 

From this it is easy to see that in case 

ii(0+) = lim < 

(including ii(0+) = — oo) we always have 



e-»0+ R{e) - 1 £ 



-00. 



In case R{0+) = 0, which means that the covariance function of Z is differ- 
entiable at the origin, define 

R"{0+) = 2 lim ^H"'^ - 
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Then we have 



C7l (0+) = Mm ^ ^ = -^y/-R"{0+). 



Hence, if i?"(0+) = — oo, which means that the covariance function of Z is 
not twice differentiable at the origin, then also C4j(0+) = — oo. Otherwise 
we are in the case where i?"(0+) is finite and Z by definition mean-square 
differentiable [l|. 

6.5 Proof of Theorem [4] 

The decompositions Et{p) = tPt{p)+Et-2{p) and Ct{p) = t'^Pt{p)+2tEt-2{p)+ 
Cffi^p) (see subsection 16.1(1 imply that ft simplifies to 

f( . _ Ctdp) Et,2{p)^ 

^'^P> Ptip) PtipY ■ 

Since C[.^{p) = Pt{p) we obtain 

We have Pt{l) = ^(t), Et{l) = ip{t), Ct{l) = tip{t) + ^(t), and 

1 



<2(i) = t / ^^^^ds+^mit)- 



Hence, ^t;2(l) = ^(t) - t^(t) and Ct-sil) = (1 + t^)^'(t) - t(p{t). By means 
of eq. (fT4l) we get 

M/2(f)-2ii;i;2(l)i5^;;2(l) 
= vI/2(t)_(^(t)_tM/(t))^(t) 



1 8 36 39 90 225 



> 0, t > 4, 



smce 

- t-'-^ + 3t"^ - Ibt-'^ < '^{t)/ip{t), t > 0. 
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Numerical inspection for t G [0,4] yields that l-2[Pt{l)]-^Et-2il)K,2i'^) 
is positive for all t > 0. Similar, but simpler argumentation holds for t < 0. 
With 



-945t" 



-11 



< '^{t)/ip{t) - (t-i - t-^ + 3t-^ - 15t-^ + lObt-^) < 0, t > 0, 



we have 

>^\t)\{e-l){t-^ 



945t 



-11\2 



2 
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2439 4935 

+ 



237825 



ti8 



+ 



tlO ' tl2 ^14 

1091475 893025 



11565 48195 
+ 



tie 



t20 



t22 



>0, 



t > 



Again, numerical inspection and similar considerations for t < yield that 
2Et.,2{lf -Ct-3{l)Pt{l)>0,teR. 

Since R"{0+) = implies that R is constant [2]J, Lemma, 1.10. 16], Lemma 
[2] yields that cov'(0+) e [-oo,0). Further, cov'(0+) = -oo if and only if Z 
is not mean-square differentiable. 
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